Numerical solutions to integral equations equivalent to differential equations with fractional time

This paper presents an approximate method of solving the fractional (in the time variable) equation which describes the processes lying between heat and wave behavior. The approximation consists in the application of a finite subspace of an infinite basis in the time variable (Galerkin method) and discretization in space variables. In the final step, a large-scale system of linear equations with a non-symmetric matrix is solved with the use of the iterative GMRES method.